Spin-dephasing anisotropy for electrons in a diffusive quasi-lD GaAs wire 
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We present a numerical study of dephasing of electron spin ensembles in a diffusive quasi-one- 
dimensional GaAs wire due to the D'yakonov-Perel' spin-dephasing mechanism. For widths of the 
wire below the spin precession length and for equal strength of Rashba and linear Dresselhaus spin- 
orbit fields a strong suppression of spin-dephasing is found. This suppression of spin-dephasing 
shows a strong dependence on the wire orientation with respect to the crystal lattice. The relevance 
for realistic cases is evaluated by studying how this effect degrades for deviating strength of Rashba 
and linear Dresselhaus fields, and with the inclusion of the cubic Dresselhaus term. 
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I. INTRODUCTION 

The ability to maximize the spin-dephasing time TjT 
of an electron spin ensemble is one of the key issues 
for developing semiconductor-based spintronic devices^'. 
However, in all III-V semiconductor materials spin en- 
sembles rapidly dephase due to the D'yakonov-Perel' 
(DP) spin-dephasing mechanism 3 -^. For the case of elec- 
tron ensembles in a heteroj unction two-dimensional elec- 
tron gas (2DEG) , two distinct contributions to DP spin- 
dephasing have to be considered: the inversion asymme- 
try of the confining potential (structural inversion asym- 
metry) and the bulk inversion asymmetry of the crystal 
lattice. The former results in an effective Rashba field 
and the latter in an effective Dresselhaus field, which in- 
cludes linear and cubic contributions 3 4 ^ 
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where x, y are the unit vectors along the [100] and [010] 
crystal directions, k x , k y are the components of the in- 
plane wave vector, and Cr_di,D3 are the spin-orbit cou- 
pling parameters. The total effective spin-orbit field B e f f 
is the vector sum of all three contributions. For 2D 
and quasi- ID electron systems, the direction and magni- 
tude of these effective spin-orbit fields can be illustrated 
as arrows on the Fermi circle. Figure [l] presents this 
for selected points in the 2D momentum space, for the 
Rashba (a) and linear Dresselhaus (b) field alone, and 
their sum (c) for the case of equal strength of Rashba 
and linear Dresselhaus field. In contrast to the individ- 
ual cases (Figure [I] (a), (b)), the magnitude of the vector 
sum shows a strong anisotropy in momentum space (Fig- 
ure [l](c)). This already suggests that spin dephasing in 
very narrow wires in which electron motion is restricted 
to the [110] direction can be strongly reduced as com- 
pared free 2DEG or such wires oriented along other crys- 
tal directions. However, it is harder to analyze whether 
such a dephasing anisotropy also occurs for wider quasi- 
lD wires, where the motion in the 2DEG plane is still 
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FIG. 1: A schematic representation of the direction and mag- 
nitude of the effective magnetic field for selected points in 
a two-dimensional fc-space, sketched for (a) the Rashba field, 
(b) the linear Dresselhaus field and for (c) the symmetric case 
of the sum of equal Rashba and linear Dresselhaus field. Both 
the magnitude and direction are depicted as arrows on the 
Fermi circle with radius &f in the {k x , fc H )-plane. 



completely random and diffusive, but where the width 
of the wire is less than length scales as the spin preces- 
sion length or the mean free path. For the latter case 
the transport regime could be named quasi-ballistic, but 
we consider the case of a large ensemble where trans- 
port along the wire is still diffusive, and where the width 
of the wire in the 2DEG plane is much wider than the 
regime of quantum confinement. Initial studies of such 
spin-dephasing anisotropy include a recent experiment 7 
on wires, and theoretical work on drifting ensembles in 
free 2DEG 8 . However, until now most emphasis was 
on work related to the spin field-effect transistor^, using 
In As-based systems or highly asymmetrical heteroj unc- 
tions, where structural i nversion asy mmetry dominates 
the spin-orbit interaction ^ ' 11 ' 12 ' 13 ' 14 ! 
We report here how the D'yakonov-Perel' spin dephas- 
ing mechanism can be strongly suppressed in diffusive 
quasi-lD electron systems based on GaAs heterojunc- 
tion material, for which Rashba and linear Dresselhaus 
spin-orbit con tributions can be indeed of comparable 
magnitude^'. The dephasing is studied for spin ensem- 
bles initially aligned perpendicular to the plane of the 
wire ([001] direction). This situation reflects the method 
of preparing and interrogating a spin population via op- 
tical pump-probe techniques 2 . Our numerical calcula- 
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tion is first performed for conditions with equal Rashba 
and linear Dresselhaus contributions and the cubic Dres- 
sclhaus term set to zero. For widths of the diffusive 
quasi- ID wires smaller than the spin precession length 
the DP spin-dephasing mechanism can be strongly sup- 
pressed and the spin-dephasing time T 2 * is considerably 
enhanced if the wire is aligned along the direction of zero 
effective spin-orbit field. Moreover, we want to point out 
that the value of our numerical tool lies in the oppor- 
tunity to study such phenomena also for more realistic 
conditions. Thus, we can study how breaking the equal- 
ity of the Rashba and linear Dresselhaus spin-orbit fields, 
or adding the cubic Dresselhaus term leads to a degrada- 
tion of the spin-dephasing anisotropy. 



II. METHOD 

We apply a Monte Carlo method^ to study the tem- 
poral evolution of the normalized spin orientation (av- 
erage spin expectation value) in an elongated quasi-lD 
wire. Our numerical tool is based on a semiclassical ap- 
proach. We use a classical description for the electron 
motion, and a quantum mechanical description of the 
dynamics of the electron spin. The wire is treated as 
a rectangular box of aspects 1 fim and 200 fj,m. The 
electron density and mobility are set to 4-10 15 m~ 2 and 
100 m 2 /Vs, which are typical values for a GaAs/AlGaAs 
heteroj unction material. All electrons are assumed to 



have the same Fermi velocity vp of 2.7T0 5 m/s. This 
is a valid approximation for fc^T, AEz.so <C Ep (with 
respect to the bottom of the conduction band), where 
AEz^so is the Zeeman splitting due to the spin-orbit 
fields alone. Electron-electron interaction and inelastic 
scattering mechanisms are neglected. 
The electron is regarded as a point particle which moves 
on a classical trajectory between scatter events on impu- 
rities (randomly determined at a rate to obtain an av- 
erage scatter time of 38 ps) yielding diffusive behavior 
in the ensemble (electron mean free path L p — vpr p — 
10 /im), and specular scattering on the edges of the wire. 
For each electron moving on such a ballistic trajectory 
we calculate the spin evolution in the effective spin-orbit 
fields quantum mechanically, and we then take the en- 
semble average on a set of electrons with random initial 
position and momentum direction. 

Within a straight ballistic segment of an electron trajec- 
tory the spin rotates around B e f f over a precession angle 
given by 
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where h is the reduced Planck's constant and t the time 
of traveling through the segment. The spin rotation oper- 
ator U for rotation over the precession angle 4> prec about 
the direction u (unit vector) of the effective magnetic field 
B e ff is obtained by (see e.gP^) 
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where a represents the vector of Pauli spin matrices 
(x, y, z components of the spin). U acting on a spin 
state Y&initiai) at the beginning of a ballistic trajec- 
tory yields the spin state \^ final) = U |* initial) at 
the end of the trajectory. Thus, we can follow the spin 
state of each electron (labeled i) , and we use this to de- 
fine a semiclassical spin vector to present its orientation, 
Si = {{S x )i, (S y )i, (S z )i) from its spin expectation values 
in x, y, z directions. 

For an electron experiencing multiple scattering events, 
the orientation of the effective magnetic field is changed 
at each scatter event. For each trajectory a rotation is 
applied. Once a scattering event takes place, the wave 
vector state is updated, i? e // is recalculated based on 
the new wave vector, resulting in a new rotation oper- 
ator, and the evolution of the spin state will carry on. 
For spin ensembles, the randomization will bring a re- 
duction of the normalized spin orientation (average spin 

expectation value) (S r ) = ,X "" X 



scmblc. (S r ) is obtained by averaging over an ensemble 
of N = 1000 spins, independent of their positions within 
the system, and calculated as a function of time. We 
choose to study the spin coherence in the ensemble here 
as (S r ). The advantage is that (S r ) gives the magnitude 
of the residual spin orientation in the direction that is 
maximum (automatically evaluating the envelope in case 
the ensemble average shows precession). However, in the 
present study without externally applied fields, there was 
no development of average spin orientation in the x and 
y directions, and the decay of (S r ) always equaled the de- 



i)/N 



The 



cay of ensemble average (S z ) 

spin-dephasing time T 2 * °f t ne ensemble is defined as the 
decay time over which (S r ) reduces to 1/e of its initial 
value. Note, however, that decay traces of (S r ) were not 
always mono exponential in our simulations. 



(Ei=i Si)/N for the en- 
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III. RESULTS AND DISCUSSION 

The temporal evolution of the normalized spin orien- 
tation (S r ) is studied for out-plane initial spin states ori- 
ented along the [001] direction. The effective spin-orbit 
field resulting from equal magnitudes of the Rashba and 
linear Dresselhaus fields is always parallel to the [110] 
axis (C R = C D1 = -1.57-10~ 8 Tmpj). These values of 
the spin-orbit parameters give rise to an average spin pre- 
cession length of about 3 /im. The width of our wire of 
1 fim is chosen to be smaller than this length scale. The 
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FIG. 2: The spin-dephasing time T% of a spin ensemble is 
plotted as a function of the wire orientation with respect to 
the [100] lattice direction. The ensemble is initially oriented 
along the [001] direction. T£ is strongly enhanced for a quasi- 
1D wire oriented in the [110] direction (black) as compared to 
an ensemble in a 2D system (gray). The arrows (in the top 
insets, horizontal axes span ±5°) for data at 45° and 225° 
indicate that here T 2 * is larger than could be calculated (Cr 
= Cm = -1.57-10 -8 Tm, Goz = 0). Inset: Ensemble spin 
expectation value (S r ) as a function of time for different wire 
orientations. 

inset of Figure [2] shows the temporal evolution for five 
different orientations of the quasi-lD wire with respect 
to the [100] direction. A distinct anisotropy of the spin 
dephasing times is observed. The peak value of T 2 is 
reached when the wire is oriented exactly along the [110] 
direction. It is found to be in excess of 10 4 ps, but the 
exact value could not be calculated in a reasonable com- 
putation time. Yet, a deviation of only 5° in the wire 
orientation with respect to the [110] direction leads to a 
reduction of T 2 * °f more than an order of magnitude to- 
wards a value of about 200 ps. For angles close to [110] ± 
15° the spin-dephasing time drops already to 43 ps and 
it reaches a minimum of 9 ps for wires oriented along the 
[110] and [110] direction. A detailed summary of these 
findings is presented in Figure[2] where the extracted spin- 
dephasing times are plotted as a function of the wire ori- 
entation. This plot reveals the strong anisotropy of T 2 
with respect to the crystal axes. 

The anisotropy of T 2 is directly connected to the motion 



of single electrons within the ensemble. Specular edge 
scattering implies that electrons with a solely transverse 
momentum component to the wire orientation (travel- 
ing less than the spin precession length between scatter 
events) almost do not contribute to the spin-dephasing 
because of motional narrowing. Only electrons with a 
strong momentum component longitudinal to the wire 
orientation are contributing to the dephasing of the spin 
ensemble. This results in the strong enhancement of Tj* 
for wires in the [110] direction. 

Spin-dephasing times for non-confined spin ensembles 
(wire width taken much wider than the mean free path 
and spin precession length) arc also calculated (gray, 
Figure [2]). In contrast to the quasi-lD wire case, no 
spin-dephasing anisotropy is found. We investigated the 
crossover from 2D to quasi ID behavior for a wide range 
of values for the spin precession length and mean free 
path (with respect to the wire width) , and found that the 
spin precession length is the crucial length scale which is 
governing this crossover. 

Next, we discuss whether this distinct spin-dephasing 
anisotropy in quasi- ID wires (with a strong enhancement 
of T 2 for wires in the [110] direction as a finger print) can 
be maintained under more realistic circumstances. First 
we study the influence of adding the cubic Dresselhaus 
term on the enhancement of T 2 * . Secondly, the influence 
of deviating strength of Rashba and linear Dresselhaus 
fields will be discussed. To avoid the difficulty that T 2 
cannot be calculated for 45° with respect to the [100] di- 
rection within reasonable calculation time, we take 43° 
as a test case, as this already shows a very strong T| en- 
hancement, while T 2 * is limited to nanoseconds. Hence, 
the calculated spin-dephasing times in the following part 
can be seen as lower bounds. T 2 * for the exact [110] di- 
rection is expected to be distinctively higher. 
The influence of the additional cubic Dresselhaus term 
on the spin dephasing time in the quasi- ID wire is pre- 
sented in Figure [3] (a). Again, the 2D case is plotted as 
a reference (grayh Without the cubic Dresselhaus term 
the calculation results in a T 2 * of 1.3 ns. However, the 
spin-dephasing time is decreasing rapidly with increas- 
ing Dresselhaus parameter. For Cd3 = -7-10 -25 Tm 3 , 
T| is already reduced to 78 ps. To estimate whether the 
spin-dephasing anisotropy is still present in a more re- 
alistic situation, experimentally deduced parameters are 
applied for comparison with our calculations. In 4 a value 
of Cd3 — -1.18T0 -24 Tm 3 is evaluated which results in 
a Tj* of only 40 ps. A value which is less than an order 
of magnitude higher than the value calculated for the 
2D case. This points to the conclusion that the cubic 
Dresselhaus term nearly annihilates the spin-dephasing 
anisotropy. However, Cd3 depends strongly on the elec- 
tron density of the system. For samples with lo wer densi- 
ties Cd3 is orders of magnitude smaller^ 17 * 18 ! . Consider- 
ing those values, it turns out, that there is a much weaker 
decay of the peak value of the spin-dephasing time, even 
when the cubic Dresselhaus term is included. Finally, T 2 
is investigated for deviating Rashba and linear Dressel- 
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FIG. 3: Data shows how the TJ enhancement of Figure [5] 
reduces when the cubic Dresselhaus term is added and when 
the symmetry between Cr and Cm is lifted. The initial spin 
state is chosen to be along [001] direction and the quasi-lD 
wire is set at 43° with respect to the [100] direction. The case 
of 43° wire orientation (rather than 45°) avoids the need to 
deal in calculations with extremely long TJ for the symmetric 
case of Cr = Cm and Cd3 ~ 0, while still clearly showing the 
ID T2 enhancement, (a) The effect of the cubic Dresselhaus 
term on TJ is plotted for the 2D case (gray) and for the quasi- 
lD case (black) (Cr = Cm = -1.57-10" 8 Tm). (b) T 2 * is 
plotted here against the difference Cr - Cm, at the fixed 
value of Cm = -1.57T0 -8 Tm for the 2D case (gray) and for 
the quasi-lD case (black) (Cdz = 0). 

haus contributions. This dependence is summarized in 
Figure [3] (b) where T| is plotted as a function of the dif- 
ference in strength of the Rashba and linear Dresselhaus 
parameter, (Cr — Cdi), in the interval ±1-10 -8 Tm for 
C m fixed at -1.57-1CT 8 Tm. At C R - C D1 = this 
results in the previously calculated T| of around 1.3 ns. 
With either increasing |Cr| or increasing |Cdi|, T£ is de- 



caying equally fast. Cr and Cdi taken fronP result in a 
difference Cr — Cdi of 0.4-10~ 8 Tm. For this value our 
calculated T 2 * is already considerably reduced to about 
110 ps. However, this spin dephasing time is still an or- 
der of magnitude higher than the one resulting from the 
2D case and in addition, Cr can be tuned with a gate or 
heterostructure design to equalize it to the linear Dres- 
selhaus field. Therefore, in summary, it can be stated 
that the spin-dephasing anisotropy which is very distinct 
for Cr and Cdi exactly equal can still prevail under less 
ideal conditions. 



IV. CONCLUSION 

A useful numerical tool is developed for studying spin- 
dephasing in device structures due to the D'yakonov- 
Perel' spin-dephasing mechanism. The Rashba, linear 
and cubic Dresselhaus contributions can be taken into 
account. With this tool is was demonstrated that quasi- 
lD wires (narrower than the spin precession length, but 
with diffusive 2D motion for the electron ensemble) can 
show very clear signatures of spin-dephasing anisotropy, 
with a strong suppression of spin dephasing for wires in 
the [110] crystal direction. 
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